Incipient failure in sandpile models 
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Elastoplastic and constitutive equation theories are two approaches based on very different as- 
sumptions for creating a continuum theory for the stress distributions in a static sandpiie. Both 
models produce the same surprising prediction that in a two dimensional granular pile constructed at 
its angle of repose, the outside wedge will be on the verge of failure. We show how these predictions 
can be tested experimentally. 
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1. Introduction 

Understanding the stress distributions in static piles is 
a basic and open question in the field of granular materi- 
als. If one tries to construct a continuum coarse-grained 
description of such systems, the difficulty is that, unlike 
for elastic materials, there is no obvious relation connect- 
ing the stresses to any strain fields. In two dimensional 
systems, there are three variables, the components of the 
symmetric stress tensor, but only two force balance equa- 
tions in the two orthogonal directions. If the three com- 
ponents of the stress tensor can no longer be expressed in 
terms of the two components of a displacement field, one 
needs one additional equation to obtain the stresses. In 
general three dimensional systems, the situation is even 
worse: one has six components of the stress tensor instead 
of three components of a displacement field, so that there 
are three 'missing' equations. 

One approach for understanding the general structure 
of the stress distributions in granular systems is in terms 
of a grain by grain statistical method However, one 
would hope that, on length scales much larger than the 
individual grain size, it should be possible to construct 
a continuum theory for a static pile of granular mat- 
ter (a "sandpile"). Several very different approaches for 
constructing such a theory have been tried. Two par- 
ticularly interesting recent models are the elastoplastic 
theories epitomized in the work of Cantelaube and God- 
dard j| and the constitutive equation models introduced 
by Bouchaud et al. ||, and subsequently generalized by 
Wittmer et al. ||. Despite their very different assump- 
tions, both approaches obtain very similar results for the 
stress distribution in a two-dimensional triangular sand- 
pile obtained by pouring grains from a point source. In 
particular, they predict that for such a sandpile there is 
a wide wedge of material near the surface of the pile that 
is on the threshold of failure, in that it is unstable to the 
application of further shear stress. 

On the face of it, the result of a large wedge at the 
instability threshold predicted by both approaches is sur- 



prising since there is a great deal of experience with the 
flow caused by tilting a pile at its angle of repose. In such 
avalanches the motion is confined only to grains in a thin 
boundary layer near the surface H |6). From such ob- 
servations one might have naively expected that only the 
surface of the sandpile, and not an entire wedge, would 
be unstable. 

These are dramatic predictions that are inescapable 
consequences of the assumptions of these two types of 
theories. It is important to see if experiments can be 
designed to detect the unstable regions. The purpose of 
the present paper is to explore the general consequences 
of having a wedge of plastic material near the surface of 
the pile, and show what experiments can be performed 
to probe its existence. 

In the constitutive equation models || @] it is sug- 
gested that the components of the stress tensor are actu- 
ally related to one another. As mentioned above, in two 
dimensions one needs one such "constitutive equation" 
interrelating the stresses; in general three-dimensional 
systems one needs three constitutive equations, which 
can be reduced in systems with special symmetry. Al- 
though the constitutive equations in general depend on 
how the system is assembled, the hypothesis that these 
equations are local implies that they are determined in 
any part of the pile when that part is constructed, and 
are not affected by subsequent changes in the loading or 
by rearrangements in other regions of the pile. Focusing 
on a two dimensional sandpile that is built by pouring 
grains from a point source, through scaling and sym- 
metry arguments, Wittmer et al. [Q are led to consider 
constitutive equations that impose a linear relationship 
among the stress components. One is ultimately left with 
a one-parameter family of models and solutions thereof. 

For a two dimensional sandpile at its angle of repose, 
all the models in the one-parameter family obtain stresses 
that are linear functions of the spatial coordinates, and 
predict that outside a symmetric triangular "inner re- 
gion" , the entire pile is on the verge of failure. This is in 
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the sense that, if <j> is the angle of repose, and r and z are 
the coordinate directions as shown in Fig. ||, the stresses 
in the outer region satisfy the condition 

*"> 2 + 4g " = Bin 2 ^ (1) 

(p rr + a zz y 

i.e. the Mohr-Coulomb criterion is saturated. This is 
equivalent to the statement that, at any point in the 
outer region, one can find a suitably oriented plane for 
which the ratio of the shear stress to the normal stress is 
tan <fi, which is the maximum possible value before plastic 
flow takes place. 




FIG. 1. Two dimensional sandpile at its angle of repose. 
The point (r 0j zo) will be used later in the paper. 

The stress distribution obtained by Cantelaube and 
Goddard || is similar to that obtained by Refs. ||] 
and M, although these authors start from a completely 
different point of view. They assume that a granular 
material can be treated as an elastic medium, with the 
added feature that when the ratio of the shear strain to 
the normal strain goes beyond a critical value, the cor- 
responding ratio of stresses saturates at tan0. Thus a 
sandpile can be divided into elastic and plastic regions, 
with the stresses expressible in terms of displacement 
fields in the elastic regions, and satisfying the Mohr- 
Coulomb yield criterion in the plastic regions. For a 
two-dimensional system, one can eliminate the strains 
from the problem and obtain expressions entirely in terms 
of stresses. (One still needs to specify the boundaries 
between the elastic and plastic regions in order to de- 
termine the stresses.) If one assumes a scaling form 
(Jij = zfij(r/z) for the stresses, one finds that the region 
adjacent to the sandpile surface must be plastic. Thus 
the elastoplastic models [|| and the ones proposed by 
Bouchaud et al. ||] and Wittmer et al. M are in agree- 
ment in predicting a wide wedge of material near the 



surface of a sandpile at its angle of repose that is on the 
threshold of failure. However a surprising, and we think 
counterintuitive, result of the elastoplastic approach (|] 
is that even for a triangular sandpile with arbitrary slope 
including one well below its angle of repose, the plastic 
region must persist if the stresses satisfy the scaling form 
above. 

In the rest of this paper, we find the destabilizing 
perturbations that can experimentally probe the plas- 
tic region in these two types of theories. Because their 
approaches are different, they must be considered sepa- 
rately. In the scenario of Refs. || and Q, an appropri- 
ately chosen infinitesimal perturbation anywhere in the 
outer region of the sandpile destabilizes it. We show this 
by demonstrating that, when the extra perturbation is 
included, there is no solution to the equations for the 
stresses that satisfies the Mohr-Coulomb yield criterion 
everywhere in the pile. The effects of such a perturba- 
tion necessarily extend to the outer surface of the pile. 
For elastoplastic theories, an infinitesimal force applied 
in an appropriate direction to a section of the bottom 
of the plastic part of the pile is sufficient to destabilize 
it. We discuss how these results would be modified for a 
material with a non-zero Bagnold hysteresis angle. 
2. Constitutive equation models 

In the scenario of Refs. j|] and |Q, if a localized body 
force of magnitude / is applied in the positive r direction 
at some point (ro, zq) in the outer region (see Fig. [[]), the 
equations satisfied by the change a' in the stress tensor 
a are 

d r a' rr + d z a' rz = fd(r - r )S(z - z a ) 

d r a' rz + 8 z a' zz = (2) 

with the boundary conditions that the change in the nor- 
mal and shear stress on the sloping sidewalls of the sand- 
pile must be zero. As discussed in the previous section, 
the arguments suggesting a local constitutive equation 
imply that this equation is not affected by subsequent 
changes in the loading, so that the components of the 
change in the stress tensor must also satisfy the equation 

Kr = W'zz + VKz ( 3 ) 

where r\ and /x are related, so that there is one free pa- 
rameter in the equations. [Q For /i ^ 0, Eq.(||) is actually 
two different equations on both sides of the symmetry 
axis r = 0, leading to a kink in the stress tensor. We 
shall ignore this at present, since we shall later see that 
it is of no consequence. 

From Eq.^), a' zz = —d r V and a' rz = d z V, where V 
is a scalar field. Substituting in Eq.([|a), we obtain 

- r$V + fxd r d z V + d 2 z V = S(r - r )S(z - z )f. (4) 

The homogenous version of this equation is an asymmet- 
ric wave equation Q] with wave velocities dr/dz = c\ and 
— | C2 1 - Eq.([|) has the solution 
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V(r,z) = \\O(c 1 {z-z )-r + r ) 

- &(\ca\(zo - z) -r + r ) 
with A = //(ci + | C2 1 ) - Therefore we obtain 
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a' zz = A <5(ci(z - z ) - r + r Q ) - <5(|c 2 |(z -z)-r + r )] 
<r' zr — A C\5{cx{z — zo) - r + r ) + 

|c2|<5(|c 2 |(2 - z) - r + r ) 

a 'rr = Wzz + Va'rz- ( 6 ) 

All components of the change in the stress tensor thus 
consist of a pulse travelling outward at a velocity c\ and 
another travelling inward with a velocity C2, as one pro- 
ceeds downward into the pile. This is as one would expect 
for a wave equation. The inward pulse is partly reflected 
and partly refracted when it reaches the symmetry axis 
r = 0, emerging as two outward pulses both with speed 

Ci. 

Since we wish to demonstrate a violation of the Mohr- 
Coulomb criterion, it is sufficient to do this for one of the 
pulses. Accordingly, we restrict our attention to the path 
of the outward propagating pulse. Along this path, we 
have 



°'zz = v'zr/ c l 



with a' zr > for / > 0. In the last of these equations, we 
have used the relation ||] r]/c\ + fi = c\ which is valid 
for all [i. When a' is small [jij, the requirement that the 
Mohr-Coulomb criterion should be satisfied is equivalent 
to the condition 

(g zz - <T r r)(P Z z - a 'rr) + ^rzKz ^ ■ 2 i ,o\ 

-— — < sin <p. (8) 

(a zz +a rr )(a' zz + o' rr ) 

One can verify that, irrespective of fi, the unper- 
turbed stresses satisfy the conditions |i) [g zz — cr rr ) = 
2<7 rz tan<^ — sin 2 4>{a zz + oy r ). Substituting in Eq.(||), 
the inequality simplifies to 



o' rz - (tan 4>)a' rr < 
which from Eq.(Q) is equivalent to 
[1 — ci tan (f>]u' zr < 0. 



(9) 



(10) 



Since a' zr > for / > 0, and c\ tan <fi < 1 in order for the 
existence of an outer region in the sandpile, we see that 
a horizontal outward directed force added anywhere in 
the outer region of the sandpile results in there being no 
possible solution to the stress equations, i.e. the sandpile 
is destabilized. 



3. Elastoplastic models 

In the scenario of Ref . , one assumes that the sand- 
pile can be divided into elastic regions, where the stresses 
are related to underlying strain fields, and plastic regions 
where Eq.(Q) is satisfied. The specific solution for a sand- 
pile at its angle of repose in Ref. || has an outermost 
plastic region in which the stresses are linear functions 
of the spatial coordinates. This is, in fact, a general con- 
sequence of assuming a scaling form Oij — zfij(r/z) for 
the stresses, valid even when the angle of the sandpile 
is below the angle of repose. Since this is not proved in 
Ref. 0, we outline the proof here. 

Choosing units in which the gravitational force per unit 
volume pg is unity, the force balance equations for the pile 
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r (T rr -\- 3 z (J rz 
d r <Jrz + d z (Tzz 



= 
= 1. 



With a scaling form for the stresses, these yield 

frr frz ~ S f rz — 
f'rz + /zz — Sfzz = 1 



(11) 



(12) 



where S = r/z is the argument of the scaling functions 
and the prime indicates differentiation with respect to 
S. In an elastic region, the fact that the stresses are re- 
lated to underlying strain fields requires the consistency 
condition 



dzdzdrr + d r d r a zz = d r d z a r 



(13) 



which with the scaling form for the stresses is equivalent 
to 



Q 2 f" _|_ f " = Q t'< 

w J rr ' J zz J r 



(14) 



Differentiating Eqs.(|12|), together with Eq. (Jl4|), we have 
three linear homogenous equations in the three second 
derivatives /y. Since the condition for the three equa- 
tions to be linearly dependent is 1 + S 2 + S 4 = 0, which 
is not satisfied for any real S, it follows that must 
all vanish, i.e. the stresses in any elastic region are lin- 
ear functions of the spatial coordinates. Since all the 
stresses must vanish on the outer surface of the sandpile, 
it follows that if the outermost region were elastic, the 
stresses would be of the form <7y = fcy-z(l — S/Sq), where 
r / z = Sq on the outer surface of the sandpile. Thus the 
ratios of the stresses to each other are independent of 
spatial coordinates in this region. From Eq.(|l|), it is then 
clear that if the Mohr-Coulomb criterion is satisfied as 
an inequality (as it must inside an elastic region), it will 
be satisfied as an inequality even on the inner boundary 
of the region, so that one cannot patch the solution to 
a plastic region. Since it is impossible |2] to obtain a 
fully elastic solution for the stresses in a sandpile, the 
outermost region must be plastic. 
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In a plastic region, Eq.(jTJ) will not be satisfied. How- 
ever, Eqs.(|l^) are still valid. For the outermost region, 
this implies that the scaling functions fij(S) are all lin- 
ear in S — So close to the surface of the pile. The pro- 
portionality constants can be obtained from Eq.([l]) with 
Eqs.(||). But since Eqs.@ with the derivative of Eq.(§) 
constitutes a system of first order differential equations in 
the scaling functions /y, one can evolve them to obtain 
the stresses everywhere in the outermost plastic region 
once they are known at the surface. Since choosing 
to be a linear function of S — So can be seen to be a valid 
solution not just close to the surface of the pile, but even 
as one proceeds into its interior, it follows that this is the 
correct solution. 

Thus we see that for any angle for the sandpile, the 
outermost region must be plastic, and the stresses are 
linear homogenous functions of (z — r/So). The extent 
of this plastic wedge is a free parameter in elastoplastic 
theories, partly constrained by the requirement of match- 
ing to an elastic region. Since the ratios of the stresses 
to each other are independent of location in the plas- 
tic zone, it follows that the yield lines, along which the 
system destabilizes under infinitesimal extra shear, have 
the same orientation everywhere in the outermost plastic 
zone. For the case when the pile is at the angle of repose, 
the yield lines are vertical and along the surface of the 
pile. If the plate supporting it at the bottom is made 
of segments, and an infinitesimal extra upward force is 
applied to a segment below the plastic region of the pile, 
the segment will move upwards. If the extra force is 
kept constant, independent of the displacement of the 
segment, the sandpile will not be able to resist the force, 
and the segment will continue to move upwards, result- 
ing in overflow at the top surface. This is in contrast to 
the response that one would obtain in the elastic regions, 
where an infinitesimal extra applied force would cause 
a slight upward motion, at which point the (deformed) 
sandpile would be able to resist the extra force. Note that 
such an argument relies on the notion that applying the 
extra force will cause the material in contact with the 
segment that is pushed to move upwards, propagating 
the disturbance upwards without attenuation. Implicit 
in this argument is the idea that stresses can be related 
to strains; since the approach of Refs. H and M does 
not consider strains at all, it is not clear that the same 
argument would apply there, which is why we obtained 
a different experimental test (in the previous section). 

Away from the angle of repose, the yield line in the 
outer plastic region in elastoplastic models tilts outwards 
from the vertical, at an angle that can be calculated in 
terms of the angle of the pile and the angle of repose. The 
above argument would apply there too, with the obvious 
modification that the extra applied force would have to 
be oriented along the yield direction. 
4. Experimental Consequences 

We now examine the effect on experiment of various 



approximations we have made: setting the Bagnold hys- 
teresis angle to zero, applying an infinitesimal force to 
a perfectly continuous medium, and considering a two- 
dimensional sandpile. 

Real granular materials do not have a well defined an- 
gle of repose. If one pours grains slowly from a point 
source, the pile slowly builds up to a maximum angle; 
adding any more grains causes an avalanche in which the 
pile suddenly collapses to an angle that is about ten per 
cent less. If one keeps adding grains beyond this, the pile 
builds up once again to its maximum angle. The Bagnold 
hysteresis angle is the difference between the two angles 
before and after the avalanche. We have treated the angle 
of repose as determining the maximum shear stress that 
any part of the pile, regardless of its past history, can 
ever sustain. One should therefore build up the pile as 
close as possible to its maximum angle before putting on 
the perturbing force. Since the pile will never be exactly 
at its maximum angle, one would always need a finite 
force to destabilize it. However, a sufficiently large force 
applied anywhere in the pile will always destabilize it, 
so that identifying a region of incipient failure becomes 
difficult. 

One can resolve this issue statistically: if the pile is 
repeatedly brought close to its maximum angle, and a 
small force is applied somewhere in the region of incipi- 
ent failure, the probability of destabilizing the pile should 
be independent of where the force is applied. One should 
be able to detect the position of the line separating the 
stable region from the outer, unstable, region by applying 
the same magnitude perturbation at different distances 
from the surface. Throughout the outer region, the per- 
turbation should be equally likely to make the pile fail. 
As soon as the position of the perturbation crosses into 
the inner region, the perturbation will become progres- 
sively less likely to cause failure as the distance from the 
separating line increases. Thus, if only a thin boundary 
layer is on the verge of flowing, the probability of desta- 
bilizing the pile will decrease rapidly when the point of 
application of the force is moved away from the surface. 

The applied force should also be sufficiently small so as 
not to destabilize a region well below the Mohr-Coulomb 
yield point; it is preferable to distribute the force over 
several grains, since a highly localized force will produce 
much greater stresses in its vicinity. Finally, if one con- 
structs a two-dimensional pile by pouring grains between 
parallel plates, the friction between the plates and the 
grains should be sufficiently small so as to avoid shear 
forces at the boundary plates stabilizing the pile. 

These issues are much less important for elastoplastic 
theories, where the extra force is applied across an ex- 
tended segment of the bottom plate, and the experiment 
can be performed even away from the maximum angle of 
repose. Of course, one would need to know the angle of 
the pile quite accurately in order to apply the perturbing 
force in the correct direction. 
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In conclusion, we have shown that the existence of 
a wide failure zone adjoining the surface of a two- 
dimensional sandpile built to the angle of repose is 
amenable to experimental verification. The presence of 
such a zone is perhaps surprising. That it should be pre- 
dicted by two very different theories is remarkable and 
deserves a careful experimental study. 
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